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ABSTRACT 

We speculate that the dark halos of dwarf galaxies and low surface brightness galaxies soften their 
central cusps by the decay of a fraction of cold dark matter (CDM) particles to a stable particle with a 
recoiling velocity of a few tens km s~ ^ , after they have driven the formation of galactic halos. This process, 
however, does not necessarily produce a significant reduction of the central dark matter density of satellite 
dwarf spheroidals like Draco or Fornax. It is shown that the recovered distribution of concentration 
parameters c for the initial (before decay) Navarro- Frenk- White halos, is in good agreement with CDM 
predictions. Other interesting potentials of unstable CDM are highlighted. 

Subject headings: galaxies: halos — galaxies: kinematics and dynamics — galaxies: spiral — dark 
matter 



1. INTRODUCTION 

The hierarchical clustering cold dark matter (CDM) sce- 
nario for structure formation has been very successful in 
reproducing most properties of the local and high-redshift 
Universe. However, some few perceived problems of such 
models, particularly at small scales, remain open. Nu- 
merical studies of coUisionless CDM with different initial 
density fluctuation spectra predict dark halos with steep 
cusps of the form p{r) ^ Po{r/rs)~°'[\ + r/rs)°'~^, with 
a = 1-1.5, after Navarro, Frenk & White (1996) (NFW 
profile, hereafter), whereas most of the rotation curves of 
dwarf galaxies and low surface brightness (LSB) galaxies 
suggest that their halos have constant density cores (see de 
Blok, Bosma & McGaugh 2003, and references therein). A 
possibly related problem is the apparent overprediction of 
the number of satellites with circular velocities Vc ~ 10-30 
km s~^ within the virialized dark halos of the Milky Way 
and M31 (Klypin et al. 1999; Moore et al. 1999). 

Many astrophysical mechanisms have been suggested to 
explain these discrepancies. If winds driven by vigorous 
star formation, bars in forming galaxies and the merger 
of black holes prove to be inefhcient at altering a primor- 
dial dark matter cusp (see Gnedin & Zhao 2002, Sellwood 
2003 and Read & Gilmore 2003, respectively), one should 
invoke modifications of the standard CDM scenario. A 
reduction of the small-scale power by substituting warm 
dark matter for CDM or by appealing to a model of infla- 
tion with broken scale invariance is strongly discouraged 
by the recent WMAP observations (Bennett et al. 2003; 
Spergel et al. 2003). 

Other alternatives involve altering the fundamental na- 
ture of the dark matter itself (Spergel & Steinhardt 2000; 
Kaplinghat, Knox & Turner 2000; Gen 2001a; Giraud 
2001; Arbey et al. 2003, and references therein). All 
these models have free parameters tunable to fit the ob- 
servations. For instance, the self-interacting dark mat- 
ter (SIDM) scenario with a general velocity dependent 
cross-section per unit mass ctd/tod = ctq {vioo/v)°' , where 
"^loo = 100 km s~^, satisfied all the astrophysical require- 
ments for a very narrow range of parameters (Tq = 0.5- 
1 cm^ g^^ and a = 0.5-1 (Hennewai & Ostriker 2002; 



Colin et al. 2002, and references therein). This huge cross- 
section has made SIDM unappealing to most WIMP and 
axion theorists. 

There are some empirical evidence against SIDM. The 
central densities of dark matter in mass models with dark 
matter cores for the elliptical galaxy NGC 4636 are higher 
than expected in the SIDM scenario with a < 1 (Loewen- 
stein & Mushotzky 2002). In addition, the halo eccen- 
tricity derived for NGC 720 is in good agreement with 
coUisionless dark matter (Buote et al. 2002). Moreover, a 
significant reduction of the central condensation of satellite 
subhalos by dark matter collisions, assisted by tidal strip- 
ping, would be difficult to reconcile with the large stel- 
lar velocity dispersions observed in the dwarf spheroidal 
galaxies (Stoehr et al. 2002; D'Onghia & Burkert 2003). 

The problems discussed above are also faced by other 
models whose mechanisms to produce a core are based on 
interactions, either elastic or inelastic, between dark mat- 
ter particles. The case of annihilating dark matter is an ex- 
ample. This Letter is aimed at investigating a new variant 
of unstable dark matter, which corresponds to a modifica- 
tion of the proposal by Gen (2001a), to find out whether 
it is possible to overcome the difficulties associated with 
previous scenarios but retaining the desired phenomena. 

2. UNSTABLE COLD DARK MATTER AND IMPLICATIONS 

Although fits to the rotation curves of LSB galaxies fa- 
vor halos with cores, a fraction of them is acceptably repro- 
duced by divergent NFW profiles as long as they are much 
less concentrated than typical halos in standard CDM. 
Gen (2001a) reconsiders the decaying dark matter sce- 
nario (e.g., Davis et al. 1981; Turner, Steigman & Krauss 
1984; Turner 1985; Dicus & Teplitz 1986; Dekel & Piran 
1987; Sciama 1990), and suggests that the concentration 
of dark matter halos lowers if one-half of the CDM par- 
ticles decay to relativistic (but nonradiative) particles by 
redshift z = 0. The decaying CDM scenario (DCDM, here- 
after) would provide also an explanation for the formation 
of dwarf spheroidal galaxies with low velocity dispersions 
(Gen 2001b). 

Unfortunately, there are some concerns associated with 
the DCDM model. First, not all the galaxies can be fitted 
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with a NFW profile varying the concentration parameter, 
especially dwarf galaxies (e.g., Blais-Ouellette, Amram & 
Carignan 2001). Secondly, there are empirical evidence 
to believe that the density of the dark matter halo of the 
Milky Way at z = 2 agrees with the present day halo 
within 50 kpc (e.g., Zhao 2002, and references therein). 
The spatial distribution and dynamics of globular clusters 
give support to this fact. On the other hand, a reduction 
in the central density of clusters of galaxies and ellipti- 
cals by a factor 2*, as suggested in the DCDM scenario, 
might be in disagreement with observations (e.g., Arabad- 
jis et al. 2002; Loewenstein & Mushotzky 2002), and could 
spoil the good agreement found in standard CDM simula- 
tions of halo substructure (Stoehr et al. 2002). However, 
we point out here that the excessive reduction of central 
dark matter densities in the halos of galaxies and clusters 
of galaxies can be alleviated if a fraction / of the CDM 
particles decays to relativistic and nonradiative light par- 
ticles plus a stable massive particle with a recoiling ve- 
locity in the center of mass of the parent particle of, say 
30-60 km s~^. If the unstable parent particles are de- 
noted by D and the stable massive particles by D', the 
decay is given by _D ^ _D' + light particles. We will re- 
fer to this model as the unstable CDM model (UCDM). 
Initially, dark halos are made up of a 'primordial' mixture 
of particles D and D' . Except for dark matter particles 
in the outer halo where the escape velocities are smaller 
than the recoiling velocity, I'rc, the particles D' created 
by decay, will remain bound to the galactic halos. Only 
the light particles that are relativistic will escape from the 
host halo. Therefore the total mass in a Milky Way halo 
within 50 kpc would remain approximately constant (see 
Eq. 2-192 of Binney & Tremaine 1987), consistent with ob- 
servations. If mo and mo' are the particle mass of D and 
D' , respectively, a recoiling velocity of a few tens km s~^ 
implies (to^ —mD')/'mD < 1- 2 x 10~'*. The consequences 
of UCDM will be discussed in the remainder of this Letter. 

2.1. Formation of Cores in Low-Mass Galaxies 

This section is sought to outline how the profile of dark 
halos would evolve when CDM particles decay and to eluci- 
date whether a central density core is expected. For sim- 
plicity we will assume pure CDM halos and ignore the 
contribution of the baryonic component. This is a good 
approximation for galaxies dominated by the dark matter. 
The implications of the inclusion of the baryonic compo- 
nent is discussed in §2.2. As long as a certain dark halo 
formed at high redshift and the mean lifetime of dark par- 
ticles is long enough, the initial halo made up of a mixture 
of particles D and D', should follow a NFW profile with 
an inner slope a « 1-1.5 and a concentration, Ci, as pre- 
dicted in cosmological simulations (Navarro et al. 1996, 
1997; Power et al. 2003). If fo is the 'primordial' fraction 
of particles D and foD'it) is the fraction of particles of 
type D that have decayed at time t, then the fraction of 
decayed particles is f{t) = fofDD'it)- Note that we do 
not need to employ a fine-tuning between the mean life- 
time of D particles, t, and the present Hubble time tn; the 
case T much shorter than <h just implies that all particles 
D have already decayed at the present Hubble time, i.e. 
fitu) = /d. 

The kinetic energy of the daughter particles D' just af- 



ter decay depends on the angle between the velocity of the 
parent particle and the recoiling velocity; some of them 
can even lose energy if, for instance, the recoiling veloc- 
ity is opposite to the direction of the orbital velocity of 
the parent particle (except when fj-c ^ o"t, where ct de- 
notes the typical one-dimensional velocity dispersion of 
halo particles). However, the kinetic energy of the daugh- 
ter particle, just after decay, averaged over all the possible 
angles of decay of a parent particle with velocity vd is 
(E) = mD'W^/2 + rriD'^^'pc/^- Therefore, the net increase 
of energy is mu'V^^./2, which is independent, on average, 
of the initial velocity v^. 

The dark matter decay will result in a gradual expan- 
sion of the halos. The typical variation in the size of halos 
comes from virial arguments as follows. A halo with a 
virial radius i?vir,i will expand after the decay of dark par- 
ticles to a virial radius i?vir,f of order of Rvii.f/Rvir,i ~ 
(1 — /.T^/3)~^, where x = frc/ct- For x small, as occurs 
in clusters of galaxies (typically x « 0.04) and giant el- 
lipticals, the effect is negligible, i?vir,f ~ -Rvir,i[l + 0{x'^)], 
even if / = 1. Accordingly, the profiles of massive galaxies 
and galaxy clusters should not deviate significantly from 
the predictions of standard coUisionless dark matter in the 
UCDM hypothesis. However, for objects with low at, as 
occurs in dwarf galaxies, the inclusion of dark matter de- 
cay may have important consequences. 

The detailed evolution of the density profile of a cer- 
tain halo will depend on the distribution function of its 
particles. The least favorable case for core formation oc- 
curs when the halo particles are assumed to be initially on 
circular orbits because this situation minimizes the kinetic 
energy in the inner halo. By contrast, in a halo dominated 
by radial orbits, the circularization of orbits by the decay 
could partly remove the cusp. In order to quantify the ef- 
fect of different internal kinematics of the CDM particles, 
we explore first the case in which most of the dark parti- 
cles are on nearly circular orbits, followed by an isotropic 
model. 

To illustrate how cores are produced, let us assume that 
the inner circular velocity Vc{R,t) of the initial configu- 
ration of the halo, when almost no particle has decayed 
yet, is given approximately by the power law Vc{R,0) = 

vo {R/ Rq)^^^'^ , with /3 > and vq the circular velocity at 
a galactocentric radius Rq. In the NFW profile, the ex- 
ponent 1/(2/3) « 1/2, whereas for a more shallow rotation 
curve, as the rigid-solid rotation, 1/(2/3) ~ 1. Suppose 
that between t = and t ~ dt a small fraction of particles 
decays (/ <C 1). The daughter particles 'gain' an amount 
of energy of AE « as said before. If we define 

the mean radius for non-circular orbits as the radius of a 
circular orbit with the same energy, it is straightforward to 
show that the mean radius of the daughter particle, Rw, 
is: 



where Rj) is the mean radius of its parent particle. A mea- 
sure of the expansion of the halo, when the initial orbits 
are approximately circular, is given by the relative varia- 
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tion in the radius of the particles: 




(2) 

which goes as cx R-q^^^ in the particular case that v^c ^ vq 
and i?D ^ Rq, keeping only the leading term. This re- 
sult suggests that the rotation curve becomes more shal- 
low, once the fraction of mass that decays is not longer 
negligible at say t = St ^ . Simultaneously, the halo will 
undergo an adiabatic expansion as a response to the new 
mass configuration, decreasing the velocity dispersion of 
the halo particles. This expansion also contributes to the 
formation of a shallow profile. Inserting the new values of 
P and vo a.t t = dt in Eq. (2), we see that 1/(2/3) always 
increases in time, concluding that the formation of a core 
should be a natural consequence in this model. For a more 
realistic, isotropic orbital distribution, the formation of a 
core should be even more favorable. 

The procedure outlined above allows us to calculate nu- 
merically the evolution of a dark halo starting with a NFW 
profile {a = 1) and for different values of v^-c (see Fig. 1). 
For the isotropic models, wc have taken advantage of the 
analytical fit to the distribution function given by Widrow 
(2002). The initial halo density profile is defined by the 
concentration parameter c and the rotation velocity V200 
at radius R200, that is the radius inside which the average 
overdensity is 200 times the critical density of the Uni- 
verse. The resulting rotation curves should be regarded 
as lower limits since the orbits of the daughter particles 
were assumed to be circular. However, the maximum in- 
crement of the rotation velocity by including an isotropic 
orbital distribution for these particles can be estimated 
from the differential energy distribution (e.g., Binney & 
Tremaine 1987) and it is marked with error bars in Fig. 1. 

In order to quantify the relevance of the UCDM model, 
we have recovered the concentration parameter before 
CDM decay, Cj, for three extreme galaxies: IC 2574, NGC 
3274 and Draco. IC 2574 is the prototype dwarf LSB 
with a conspicuous core (e.g., Blais-Ouellette et al. 2002) 
whereas NGC 3274 is a LSB galaxy with one of the small- 
est core-radius, which can be fitted with a NFW profile 
(de Blok & Bosma 2002). On the other hand, Draco is a 
dwarf spheroidal (dSph) and presents central dark matter 
densities (within a radius of ^ 1 kpc) as high as those 
predicted by standard coUisionless CDM models (Kleyna 
et al. 2001; Lokas 2002; Stoehr et al. 2002; Hayashi et 
al. 2003). Therefore, these galaxies should bracket the con- 
centration parameters for low-mass galaxies. If the initial 
conc;c!ntration parameters for these extreme galaxies, c^, 
lie within la scatter of the predictions of standard CDM 
simulations, it would imply that our model is fully consis- 
tent with the observed distribution of concentrations. For 
the reference values / = 0.5 and fj-c = 40 km s^^, the ob- 
served rotation curves are reproduced well in the isotropic 
model, starting with a NFW halo with c, « 13, 26, 20 and 
V2oo,i » 60,55,60 km s'^ for IC 2574, NGC 3274 and 
Draco, respectively. Although cosmological N-body simu- 
lations must be done for a fully description of the profile 

1 In the DCDM scenario proposed by Cen (2001a), the byproducts of 
implies that the adiabatic evolution of the halo is self-similar and he 
brealfs down the self-similarity in the evolution. 



of UCDM halos, these results are promising if one consid- 
ers that expected values of q in ACDM cosmology have 
(T(log Ci) = 0.18 and a 2a range from 5 to 25, with an aver- 
age of around 10 to 15 (NFW; Eke et al. 2001). Since the 
effects of UCDM are smaller for more massive halos, the 
present mechanism for core formation is fully consistent 
with the fact that the best-fitting models of high surface 
brightness galaxies favor cores with small sizes (Jimenez, 
Verde & Oh 2003). 

2.2. The Case of Dwarf Spheroidals 

Baryonic effects on the dark matter distribution could 
affect our estimate of Cj. In particular, if the progeni- 
tors of dSph were much more massive than presently, the 
enhancement of the velocity dispersion of CDM particles 
by adiabatic compression would have effectively reduced 
the evacuation of dark matter particles. In the following 
we will analyze how massive the dSph progenitors should 
have been to suppress significantly the effect of CDM de- 
cay. For simplicity, the baryons are assumed to dominate 
the total mass (dark matter plus baryons) within the op- 
tical radius of the galaxy during the phase of CDM decay. 
The corresponding density profile in this phase is taken to 
be p ^ R^^, with ^ w 1.8- 2, for consistency with other 
studies (e.g., Hennawi & Ostriker 2002). The flat circular 
velocity associated with this configuration will be denoted 
by Vc- Under the approximation that the CDM particles 
are on circular orbits, a CDM particle initially at a or- 
bital radius R will shrink to a radius R' = {M/M')R, 
where M is the initial mass within R and M' is the final 
mass within R' , during the dissipative collapse of baryons. 
After decay, the daughter particle will acquire a mean cir- 
cular radius R{M/M') exp (^v^^/2v'^) . Later on, most of 
the baryonic mass will be ejected probably by supernova 
explosions or stripped by ram pressure. This mass loss 
will expand the orbit of the daughter particle to a final 
radius of ^ i?exp (ti^(,/2i;^), still assuming that the pro- 
cess evolves adiabatically. The net effect is a self-similar 
expansion R XR, with A = exp (^v^^/2v^) and, there- 
fore, the initial inner velocity profile v{R) oc i?,^/^, associ- 
ated with a NFW profile pu cx R"^, is transformed into 
v'{R) = vv{R) with = 1- f + //A2. We will employ 
a variation oi y k. 0.8, which is compatible with observa- 
tional uncertainties (> 25%) in the mass determination of 
dSph (e.g., Kleyna et al. 2001; Lokas 2002), and the fact 
that, in our model, the mass loss of CDM particles induced 
by Galactic tidal forces would be significantly suppressed if 
the progenitor was more compact in the past. For / = 0.5, 
Wic = 40 km s^^ and v = 0.8, a value Vc ^ 35 km s~^ is 
required. In terms of the baryonic mass, this implies a 
mass of gas > 2.5 x 10^ Mq for the progenitor of the dSph 
Fornax with a typical radius of ~ 0.8 kpc, and > 0.4 x 10^ 
M0 for the progenitor of Draco with radial size of ^ 0.15 
kpc. As we have assumed circular orbits, these estimates 
are clearly lower limits. For progenitors with these masses, 
the decay of UCDM could not have observable effects in 
their halos. The inferred masses are large, especially for 
Draco, but are not unreasonable; they are in accordance 
with different lines of research that point out masses of 

decay were assumed to be relativistic. This assumption immediately 
no core is formed. By contrast, the inclusion of a recoiling velocity 
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10^-10^" Mq for the progenitors of those dwarf sphcroidals 
with a long period of star formation hke Draco and For- 
nax (Kitayama et al. 2001; Ikuta & Arimoto 2002). In a 
much more speculative scenario, if the gas mass of dSph 
is never larger than 10^ Mq, dSph could preserve their 
initial dark hales intact under the assumption that CDM 
only decays by interactions with baryons/nucleons with a 
certain cross-section ctdn- 

3. DISCUSSION 

The possibility that dark matter decays is being a topic 
of active work. Theoretical predictions for lifetime of 
heavy sterile neutrinos and Ultra Heavy particles, as WIM- 
PZILLAs, cover a large range of values, from 10~^ to 10^^ 
Gyr. We have assumed that UCDM and its byproducts 
are weakly interacting particles. The possibility that rela- 
tivistic light particles of the decay are the source of Ultra 
High Energy Cosmic Rays or responsible; for the reioniza- 
tion of the Universe is still open. Even if Ultra Heavy 
particles were to decay to Standard Model particles, they 
are not ruled out by present observations (e.g., Ziaeepour 



2000). 

Besides the formation of density cores, other interest- 
ing phenomena associated with UCDM arise. UCDM en- 
hances substructure destruction within galactic halos, and 
the evacuation of part of the central mass of dark matter 
in dwarf ellipticals and dwarf spheroidals would lead to an 
adiabatic expansion of their globular cluster systems, and 
would partly compensate the drag of dynamical friction, 
alleviating the problem of the survival of the globular clus- 
ters (Lotz ct al. 2001, and references therein). Moreover, 
assuming that UCDM particles are self-interacting with 
the cross-section found by Hennawi & Ostriker (2002), 
o'Do/mD ~ 0.02 cm^ g^^, it is feasible to incorporate 
the growth of supermassivc black holes in our model, as 
suggested by Ostriker (2000). The decay of particles can 
contribute to feed up the central black hole if they are 
scattered into the loss-cone, but this contribution is much 
smaller than the Bondi accretion. 

I thank V. Avila-Reese and an anonymous referee for 
helpful comments. 
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Fig. 1. — Variations in the shape of the halo contribution to the rotation curve for different values of / and Vic, starting with a NFW 
profile with ci = 13 and V^oo.i = 60 km (triple-dotted dashed line) except for the NGC 3274 halo for which Ci = 26 and V^oo.t = 55 km 
s~^. Isotropic models arc labeled by ISO, while models in which all particles are initially on circular orbits are indicated by CIR. The three 
solid lines show the rotation curves for / = 0.5 and i>rc = 40 km s~^. The dashed lines correspond to Urc = 60 km and / = 0.5 (lower 
dashed curve) and / = 0.25 (the two upper dashed curves). The bars account for the maximum effect when including non-circular orbits for 
the daughter particles for the IC 2574 and NGC 3274 halos. 



